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Abstract 
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In  1974  Doctor  David  W.  Berrie  of  Aeronautical  Systems 
Division,  proposed  a radar  jamming  waveform  derived  from  fil- 
tered pseudo-random  binary  sequences.  The  running  time  of 
the  sequence  is  synchronized  to  the  pulse  repitition  rate  of 
the  victim  radar  and  defeats  the  post  detection  pulse  inte- 
grator. Berrie 's  results  were  derived  from  theory  and  ana- 
log computer  simulation.  This  study  applies  that  waveform 
to  a rada^-  receiver  model  incorporating  a DICKE-FIX  (wide- 
band limiter)  ECCM  device  and  uses  a digital  computer  simula- 
tion (RADSIM)  to  determine  radar  performance.  Berrie 's  re- 
sults were  essentially  confirmed,  and  the  DICKE-FIX  device 
was  found  to  be  ineffective  against  this  form  of  jamming. 

/i 
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I.  Introduction 


In  197^  Doctor  David  W.  Berrie  of  Aeronautical  Systems 
Division,  USAF  Systems  Command,  proposed  a radar  jamming 
waveform,  derived  from  filtered  pseudo-random  binary  se- 
quences. The  running  time  of  the  sequence  was  to  be  syn- 
chronized to  the  pulse  repetition  rate  of  the  victim  radar 
and iwas  intended  to  defeat  a post-detection  pulse  integra- 
tor. He  was  able  to  demonstrate,  using  analog  computer  sim- 
ulation  techniques,  that  the  probability  of  target  detection 

I 

was  reduced  when  this  form  of  jamming  was  used,  as  shown  in 
Figure  1.  The  analysis  was  conducted  on  a model  of  a simple 
search  radar  with  a 127KHZ  IF  bandwidth  and  no  Electronic 
Counter  Counter  Measures  (ECCM)  devices.  The  purpose  of 
this  thesis  will  be  to  adapt  an  existing  digital  computer 
simulation  (RADSIM)  to'  analyze  a search  radar  incorporating 
ECCM  devices  and  to . demonstrate  that  synchronized  pseudo- 
random rioise  jamming  (SPRNJ)  will  be  effective  against  a 
specific  ECCM  device  known  as  DICKE-FIX,  a wide-band  limiter. 

Background 

Electronic  warefare  (E/J)  has  played  a major  role  in 
military  operations  since  the  beginning  of  World  War  II, 

! 
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Probability  of  Detection  (P 


Signal  to  Noise  Ratio  (SNR) 


Figure  1.  Relative  Effect  of  Random  Noise  And  Synchronized 
Pseudo-Random  Noise  of  Detection  Probability 
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which  introduced  the  first  large  scale  use  of  electronics  in 
military  operations  (Ref  1:1-1).  It  is  intended  to  deny 
enemy  forces  the  use  of  his  electronic  equipment,  to  use  his 
own  electronic  emissions  against  him,  and  to  protect  ones 
own  electronic  equipment  against  enemy  action.  EW  must  then 
be  highly  responsive  to  the  action  of  an  enemy  and  perhaps 
more  importantly  to  his  potential  for  action  (Ref  1:1-2). 

The  need  for  a particular  capability,  however,  only  becomes 
apparent  when  the  existence  of  a threat  has  been  established 
or  postulated.  It  is  my  own  opinion  that  the  development  of 
EW  systems  is  deficient  to  some  degree  in  the  ability  to  re- 
spond. Funding  limitations  generally  inhibit  development 
from  proceeding  beyond  the  conceptual  phase  until  a firm 
threat  has  been  established.  Thus  it  is  possible  for  an 
enemy  to  deploy  a carefully  concealed  technique  to  which  we 
would  have  no  coimtermeasure . Flexibility  is  needed  in  EW 
systems  to  meet  the  potentially  diverse  spectrum  of  threats. 

A cheap  method  of  determining  the  efficacy  of  EW  sys- 
tems against  postulated,  but  unverified  threats  is  necessary, 
since  a usual  constraint  on  development  is  funding.  In  this 
way  hardware  development,  when  it  became  necessary,  could 
proceed  with  confidence  and  rapidity.  Computer  simulation 
can  provide  a cheap  and  flexible,  in-house  analysis  tool 
for  Air  Force  research  and  development  organizations. 

Analysis  of  the  Problem 

Assumptions.  Two  major  assumptions  are  made  in  the  de- 
velopment of  this  thesis,  to  simplify  the  work.  They  are 
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not  central  to  the  main  purpose  of  developing  a tool  for  the 
analysis  of  ECM/ECCM  and  radar  performance; 

1)  Doctor  Berrie's  development  of  the  performance 
of  a post-detection  integration  radar  in  the  presence  of 
synchronized  pseudo-random  noise  is  essentially  correct. 

The  experiment  which  will  be  run  to  test  the  effectiveness 
of  SPRNJ  against  a DIGKE-FIX  radar  receiver,  will  be  de- 
signed in  accordance  with  the  procedure  he  outlines. 

2)  The  RADSIM  computer  simulation  is  correct  in 
its  signal  transfer  algorithms  and  radar  modular  transforms. 
Errors  in  these  functions,  as  they  apply  to  the  computer  to 
be  used,  will  be  corrected  within  the  context  of  the  original 
design.  The  radar  model  algorithms  are  relatively  easy  to 
change  due  to  its  modular  structure,  and  in  fact  are  speci- 
fically intended  for  modification  as  a user  requires  (Ref  3: 
^-3). 

Approach  to  the  Solution.  The  computer  simulation  to 
be  used  in  this  study  was  developed  under  a Rome  Air  Devel- 
opment Center  (RADC)  contract  for  use  on  a Honeywell  6180 
Computer  using  the  General  Comprehensive  Operating  System 
(GECOS)  and  FORTRAN  IV.  It  will  be  converted  to  the  CDC 
6600/CYBER7^  Computer  using  SCOPE  J.U  and  FORTRAN  Extended. 

Upon  completion  of  the  necessary  recoding  and  opera- 
tional debugging,  a series  of  test  cases  will  be  run  to  in- 
sure that  the  operation  of  the  simulation,  duplicates  the 
same  set  of  cases  as  run  on  the  Honeywell  6180  Computer. 

Then  an  experiment  to  simulate  SPRNJ  against  a simple  search 


r 


radar  will  be  designed  and  executed  to  confirm  the  results 
of  Doctor  Berrie's  theory  as  demonstrated  on  the  analog 
computer.  An  ECCM  device  (DICKE-FIX)  will  be  added  and  its 
performance  analyzed.  A receiver  operating  characteristic 
curve  for  the  radar  will  be  developed  for  the  two  cases; 
with  and  without  the  ECCM  device. 
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II.  The  Theory  of  Pseudo-Random  Noise  Jamming 

Doctor  Berrie's  dissertation,  identified  as  ASD  TR-74- 
23  (Ref  2),  "The  Effect  of  Pseudo-Random  Noise  on  Radar 
System  Performance"  provides  a very  readable  and  comprehen- 
sive discussion  of  the  theory  behind  the  generation  of 
pseudo-random  noise  and  its  use  as  a radar  jamming  waveform. 
This  section  will  summarize  the  theory  presented  by  Doctor 
Berrie,  and  describe  some  of  his  results  as  they  relate  to 
this  investigation. 

Pseudo-Random  Binary  Sequences 

Binary  sequences  are  sequences  which  are  two  valued, 
usually  one  and  zero.  A sequence  may  or  may  not  be  periodic, 
i.e.,  after  some  number  of  values  the  sequence  begins  to  re- 
peat. A truly  random  sequence  can  never  be  periodic,  how- 
ever, the  period  of  a periodic  sequence  could  be  so  long 
that  over  reasonable  observation  intervals  it  would  appear 
to  be  random.  If  a periodic  sequence  shares  certain  proper- 
ties of  truly  random  sequences  it  is  termed  pseudo -random. 

The  required  properties  are: 

1)  The  number  of  ones  equals  the  number  of  zeros, 
plus  or  minus  one . 

2)  of  the  total  number  of  runs  (consecutive 
elements  of  the  sequence  having  the  same 
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value)  have  length  n. 

3)  The  auto  correlation  function  of  the  sequence 
C(t)  is  two  valued  and  has  a value  C(0)  when 
T=0  or  an  interger  multiple  (k)  of  the  period 
(T)  and  some  other  value,  C(t)  when  t/O  or 
kT  (Ref  4:280) . 

The  form  of  the  function  which  generates  a periodic  binary 
sequence  determines  whether  or  not  it  is  pseudo-random.  A 
function  which  produces  a maximal-length,  linear,  periodic 
binary  sequence  will  form  a pseudo-random  sequence  (Ref  2:18). 


Generation  of  Pseudo -Random  Sequences 

In  order  to  meet  the  requirements  for  generation  of  a 
linear  pseudo-random  sequence  the  expression  for  the  recur- 
rence relationship  which  describes  the  sequence  must  be 
linear  and  of  the  form: 


^n-l  (2-1) 

aj^=  the  term  of  the  sequence  (1  or  zero) 


a binary  coefficient  (1  or  zero)  of  the 
n-i^*^  term. 

^n  i~  preceeding  term  of  the  sequence 
r=  the  number  of  initial  terms  in  the 
sequence . 

The  summation  is  a modulo  two  summation  where  the  results  of 
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the  siiinraation  are  shown  in  Figure  Z.  Equations  (2-2a) 
through  (2-2d)  provide  an  example  of  a linear  sequence 


= a,,C.+aXo+aoC..+a.  C. 

S'!  j ^ ^ J L ^ 

(2-2a) 

= a^(l)+a^(0)+a2(0)-ra^(l) 

(2-2b) 

= (l)(l)+(l)(0)+(0)(0)+(l)(l) 

(2-2c) 

= 0 

(2-2d) 

satisfying  Equation  (2-1)  with  four  initial  terms.  The 
full  sequence  will  have  2^-1  terms,  in  this  case  I5,  and 


would  describe  all  possible  states  of  an  r -bit  binary  num- 
ber less  the  all  zero  state.  The  al2  zero  state  must  be 
excluded  to  prevent  the  sequence  generation  from  stopping. 
If  Equation  (2-2a)  v/ere; 


= ai^C^+a.jC^+aL^C^+a.^Ci^  (2-3) 

= a^(l)+a_(0)+a2(0)+a^(l) 


U' 


=(0)(l)+(0)(0)+(0)(0)+(0)(l) 

=0 

no  further  terms  could  ever  be  generated. 

A generating  function  for  a linear  pseudo-random  se- 
quence may  then  be  defined  as: 


G(x)=E 


n=o 


(2-M 


where  a^  is  defined  in  (2-1).  Substituting  (2-1)  into  (2-4) 
leads  to  the  following  expression  (Ref  2:18). 


G(x)=E  C.X^E  a-kX 
^ k=l 


-k 


1-  E C^X^ 
i=l 


(2-5) 


This  may  be  redefined  as  the  ratio  of  two  polynomials. 


G(x)=- 


g(x) 

f(x) 


(2-6) 


The  polynomial  f(x)  is  the  characteristic  polynomial  of  G(x) 


f(x)=l-E  C.X^ 
i=l  ^ 


(2-7) 


r 

»1+E 


C.X" 


i=l 


If  f(x)  is  a primitive  (does  not  divide  modulo  two  into  X^+1 
for  values  of  K less  than  2^-1)  and  irreducible  (cannot  be 
factored  into  two  terms  of  lower  degree)  it  will  produce  a 
maximal  length,  linear,  pseudo-random  binary  sequence. 
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ENAD  TM-7^-01  (Ref  5)  also  produced  "by  Berrie  contains  a 
tabulation  of  primitive-irreducible  polynomials  sorted  by 
the  number  of  non-zero  coefficients.  This  tabulation  will 
be  later  used  in  developing  a shift  register  to  produce 
pseudo-random  noise  in  Section  IV. 

Generation  of  Pseudo-Random  Noise  with  Shift  Registers 

Equation  (2-1)  can  be  easily  implemented  in  a shift 
register.  Figure  3 shows  an  implementation  of  Equation 
(2-1)  for  r=4.  Given  some  non- zero  initial  state  each  stage 
of  the  shift  register  is  shifted  to  the  right  upon  a clock 
pulse  and  a^^  is  shifted  into  the  left  most  stage.  The  con- 
tents of  stages  with  non-zero  are  modulo  two  s\immed  to 
form  a new  a^ . The  sequence  may  be  obtained  one-bit-at-a- 
time  from  the  a^  stage  or  r sequential  bits  of  the  sequence 
obtained  at  one  time  from  the  r registers  stages  (Ref  2:31) • 
If  the  Cj^'s  are  determined  by  the  non-zero  terms  of  a prim- 
itive irreducible  polynomial,  the  sequence  generated  will  be 
a pseudo-random  sequence  of  length  15. 

If  the  all  zero  state  could  occur,  the  content  of  each 
stage  could  be  considered  an  independent  binary  random  vari- 
able, and  the  sum  of  all  the  stages  would  also  be  a random 
variable. 

r 

Z=E  a.  (2-8) 

i=l^ 

Z would  assume  values  from  zero  to  r with  a probability  de- 
fined by  the  binomial  distribution  (Ref  2:32). 

rN_k  r-k 


(2-9) 


Figure  3*  Shift  Register  iTiplementation  of  Linear  Recurrence  Relationship 


When  p=q=2  Equation  (2-9)  becomes 

p(Z=k)=(J)(i)^  (2-10) 

Papoulis  has  shown  (Ref  6:67)  that  for  large  r (r»4),  that 
2-10  is  asymptotically  equal  to  the  Gaussian  distribution. 

P(  Z=k)«-7~^expr -2(k-r/2)  ^/r1  (2-11) 

rrr 

with  a mean  |i=r/2  and  variance  a^=r/4.  This  is  not  strictly 
true  for  the  shift  register  of  Figure  3 since  the  all  zero 
state  cannot  occur,  but  for  sufficiently  large  r this  will 
not  significantly  alter  the  result  (Ref  6:67) » and  the  dis- 
tribution remains  essentially  Gaussian. 

The  Power  Density  Spectrum  of  Pseudo-Random  Noise 

By  Papoulis  (Ref  6:3^1)  the  auto  correlation  function 
of  a random  binary  transmission,  Figure  4a,  shown  in  Figure 
4b  is  a triangle  and  its  power  spectrum  is  defined  by  (2-12) 
and  depicted  in  Figure  4c. 

S(m)=C(o)4sin^(mT/2)  (2-12) 

——3 

The  pseudo-random  binary  sequence  is  periodic,  however,  and 
its  spectrum  will  be  composed  of  discrete  spectral  lines 
whose  number  will  be  determined  by  the  sequence  length. 

This  is  shown  in  Figure  4d.  The  expression  for  the  power 
spectrum  of  the  pseudo-random  sequence  is; 

00 

S^((t))=2'TTr 6 ((|))+N+1S  sin  n-rr/N  ® 6(a3-n2Tr)1  (2-13) 

I jp- htt/N  L 

n=-oo 

N = the  number  of  bits  in  the  sequence 
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Figure  4.  Random  Binary  Transmission,  Autocorrelation 
Function  and  Power  Spectrum 
13 
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The  zeros  of  this  spectrum  will  occur  when  n/N  is  an  integer 
multiple  of  the  clock  frequency. 

It  has  been  shown  (Ref  2:33-35)  that  the  summation  of 
the  stages  of  the  register  is  equivalent  to  low  pass  filter- 
ing of  the  sequence.  This  results  in  a decrease  in  the 
width  of  the  spectrum  main  lobe  by  a factor  of  l/r  where  r 
is,  as  before,  the  number  of  stages  in  the  shift  register. 
The  spectrum  now  appears  as  in  Figure  5* 

Generation  of  PRN  Signal  for  Jamming 

The  spectrum  of  the  shift  register  pseudo-random  noise 
generator  contains  elements  which  range  in  frequency  from 
zero  to  the  clock  frequency  divided  by  r.  To  be  useful  for 
jamming  this  spectrum  must  be  translated  to  radar  frequen- 
cies. A method  of  doing  this  would  be  to  use  the  sum  of  the 
registers  to  modulate  a sinusoidal  carrier.  In  this  study, 
double  sideband  (DSB)  modulation  will  be  used  and  Equation 
(2-14)  will  be  the  expression  of  the  jamming  signal. 

J(t)  = Z (t)cos(w  t+)2()  (2-14) 

J(t)  = the  jamming  signal 
Z = the  register  sum 

= the  carrier  frequency 
^ = some  phase  angle 

If  this  signal  is  within  the  radar  bandpass  some  degree  of 
interference  will  be  seen  in  the  radar  detection  process. 
Figure  6 shows  a block  diagram  of  the  jammer,  and  the  trans- 
mitted spectrum. 


S(w) 


Unfiltered 


Filtered 


Figure  5- 


Filtered  and  Unfiltered  Pseudo-Random 
Sequence  Power  Spectrum  Envelope 


Figure  6.  Block  Diagram  of  PRN  Jammer  and  its  Spectrum 
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The  Application  of  Pseudo  Random  Noise  to  Radar  Jamming 

A block  diagram  of  a typical  search  radar  is  shovm  in 
Figure  ?•  The  transmitter  sends  out  relatively  short  R-F 
pulses  at  specific  intervals  (typically  200-400  pulses-per- 
second) , controlled  by  a stable  local  oscillator.  The 
pulses  are  reflected  from  the  target,  received  at  the  radar 
antenna  and  reduced  to  some  convenient  intermediate  fre- 
quency (IF).  They  are  then  detected  and  used  to  drive  some 
form  of  display  which  presents  the  pulses  receivec  • . a 
scale  relative  to  the  total  round-trip  time  from  transmit- 
ter to  target  and  back  to  the  receiver. 

The  detection  process  is  essentially  a binary  detec- 
tion process  with  two  hypotheses  (Ref  8:27). 

r(t)=n(t);  noise  only  (2-15) 

r( t )=n( t)+s( t) ; signal  + noise 

r(t)  = the  received  signal 
n(t)  = noise 

s(t)  = the  reflected  signal 

The  intermediate  frequency  filter  is  a matched  filter  in 

most  radars  (Ref  7 '-7^)  whose  impulse  response  is  matched  to 

the  form  of  the  trainsmitted  pulse  S(t). 

H(w)=S*(w)e“^^‘‘^o  (2-16) 

H((i3)=filter  frequency  characteristic 
S (co)=signal  spectrum  complex  conjugate 

When  only  noise  is  considered  is  the  hypothesis  and  r(t)= 

n(t).  The  signal  out  of  the  filter  (Figure  7)  is  x('t). 

x(t)=N(t)cos[(0Qt-^(t)]  (2-17) 

N(t)  and  ^(t)  are  independent  random  variables  which  are 

normally  and  uniformly  distributed  respectively.  The  sam- 
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(2-18) 


pies  of  the  noise  are; 

W(N.)= expC-N.  V2crjj^] 

with  zero  mean  and  variance  . If  the  target  signal  is 

present,  then  is  true. 

r(t)=N(t)+S(t) 

=N(  t)+A[l+cos(i)o  t] 

=(N^+A)  cosWot+N^sintOe  t 
=Xj^(  t)  cosco^jt+X^C  t ) sinwjit 

and  X2  are  statistically  independent  random  variables 
and:  W(X^)=  1 expr-(X^-A)  (2-19) 

V2nOjg 

W(Xo)=  1 exD\- -xyza^^n  (2-20) 

W(X^,X2)=_i_expC-(X^-A)+X2^]  (2-21) 

Translating  2-21  to  polar  coordinates  by  setting  Xj^=R  cos  0 

and  X2=R  sin  0 we  can  integrate  with  respect  to  0 and  normal 
ize  with  respect  to  Ojq  to  obtain  the  joint  probability  den- 
sity function  of  the  IF  filter  output. 

W(v)=ve"^^^'*‘^^^I,,(av)  (2-22) 

a=A/0j^ 

v=R/0j^ 

Io(av)=zero  order  modified  Bessel  Function. 
Since  a is  the  amplitude  of  the  signal,  s(t),  for  the  noise 
only  case,  a in  Equation  (2-22),  will  be  zero  making  the  den^ 
sity  Rayleigh.  For  a signal  present,  a will  have  some  non- 
zero value  and  the  output  will  be  a Generalized  Rayleigh 
(Rician)  distribution  (Ref  2!l77).  Figure  8 shows  a plot  of 


Equation  (2-22)  for  the  noise  only  and  a signal  plus  noise 


case. 


Probability  of  detection  is  defined  as  the  probability 
that  the  signal  plus  no  voltage  will  exceed  some  thresh- 
old voltage  Yi  probability  of  false  alarm  is  defined  as 
the  probability  that  the  noise  alone  will  exceed  the  same 
threshold.  These  probabilities  are  shown  pictorially  in 
Figure  9*  Their  expressions  are  shown  in  Equations  (2-23) 
and  (2-24) . 


00  -v®+a^ 

P^=/ ve  I ( av ) dv 

Y 2 

00  -v^/2 


P 


dv 


(2-23) 

(2-24) 


Pulse  Integration 

Some  radars  use  pulse  integration  to  improve  the  proba- 
bility of  detection.  The  pulse  integrator  shown  in  Figure 
10  adds  several  pulses  after  detection  and  prior  to  the 
thresholding  of  Equations  (2-23)  aJ^d  (2-24).  The  impulse 
response  of  the  integrator  is: 

k-1 

h(t)=  I 6(t-kT)  (2-25) 

k-o 

K=number  of  pulses  to  be  integrated 
T=the  pulse  repitition  interval 

thus  d(t)=y( t)*h( t)  and, 

K-1 

d(t):;=  E y(t)(t-kT)  (2-26) 

k=0 


If  y(t)  is  periodic: 
d(t)=Ky(t) 
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(2-27) 
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When  y(t)  is  a random  variable,  the  sum  of  the  K additions 
will  also  be  a random  variable,  normally  distributed  with 
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mean  and  variance  . 

k 

|a^=mean  of  the  input  signal 
0^=variance  of  the  input  noise 

The  distributions  for  the  envelope  of  the  noise  only  and  the 
signal  plus  noise  are  respectively  normalized  to  the  input 
noise  variance: 

P^(v)=___l exp[-(v-K75)^/2k(2-Tr/2)]  (2-30) 

V2TrK(2-rT/2) 

and 

P_(v)__l exp[-(v-K|iy<jj^)  V2KOp®/Oj^®] 

®'^"'V2Tr(2-TT/2)  0^  (2-31) 


(2-28) 

(2-29) 


The  probability  of  detection  and  false  alarm  may  be  found  by 
integrating  Equations  (2-30)  and  (2-31)  from  the  threshold 
to  infinity  as  in  Equations  (2-23)  and  (2-24)  and  shown  in 
Figure  9. 

Figure  11  shows  the  result  of  such  integration  per- 
formed by  Berrie  for  eight  pulses  integrated  at  several 

signal- to-noise  ratios.  ; 

i 

When  synchronized  pseudo-random  noise  is  present  it 

is  a periodic  signal.  This  by  Equation  (2-2?)  its  form  will  \ 

! 

be  the  same  at  the  output  of  the  integrator  as  at  the  input.  i 

3 
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If  P(v)=ve"  lo(av) 

s+n  2 


then; 


-v^+( ak) ' 

P(kv)=kv  e 2^ 
s+n  IP 


Io(f) 


(2-32) 

(2-33) 


When  P^  and  P^^  are  calculated  for  the  integrated  signal,  no 
improvement  over  the  pre-integrated  case  is  found  (Ref  2; 180), 

Results  of  Berrie*  s Experiment 

Berrie  anailyzed  the  performance  of  synchronized  pseudo- 
random noise  against  a typical  search  radar.  He  used  an 
analog  computer  model  of  the  radar  scaled  to  audio  frequen- 
cies for  the  analysis.  The  results  of  his  experiment  showed 
that  the  probability  of  detection  was  reduced  in  a post  de- 
tection integrator  when  this  form  of  jamming  was  present. 
Table  I shows  the  results  of  his  experiment,  for  thermal 
random  noise  (TRN)  and  SPRN. 


Table  I 

Probability  of  Detection  For  P„  =0.1 
at  Integrator  Input  and  Output 


SNR 


P.  Integrator  Input  P . Integrator  Output 

^ /rm^TVT  ^ 


2 

.45/. 31 

.30/. 64 

4 

.70/.  41 

.40/. 74 

6 

.75/.  52 

.50/. 80 

8 

.85/.  65 

.60/. 84 

10 

.92/.  71 

.70/. 86 

As  would  be  expected  the  improvement,  or  rather  degradation 
in  detection  probability,  is  lessened  with  increasing  SNR, 
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effect  of  synchronized-pseudo-random  noise  based  jamming  was 
summarized.  This  treatment  was  by  no  means  exhaustive.  Most 
of  the  detail  is  treated  in  ASD  TR-74-23,  and  a reader  in- 
terested in  the  derivations  should  refer  to  that  document. 


III.  The  Radar  Simulation 


One  of  the  objectives  of  this  study  is  to  provide  to 
AFIT  and  ASD  a computer  program  for  the  analysis  of  radar/ 
communications  hardware  and  waveforms.  This  objective  is 
met  by  converting  an  existing  program,  developed  by  the  Rome 
Air  Development  Center  for  use  on  a Honeywell  Information 
System  635/6I8O  computer  to  the  ASD  CDC  6600/CYBER74  com- 
puter. This  section  will  describe  the  development  of  the 
RADSIM  simulation  program,  and  explain  the  major  changes  re- 
quired to  convert  it  to  the  CDC  computer.  It  will  also 
briefly  explain  how  to  structure  a problem  and  prepare  data 
cards.  Complete  program  documentation  is  provided  in  the 
four  volumes  of  RADC  TR-76-I86,  and  in  the  source  code  as 
altered  to  rim  on  the  CDC  system.  While  the  simulator  main 
programs  execute  correctly,  some  of  the  subroutines  remain 
unverified  for  CDC  use. 

Background  of  RADSIM  Development 

The  RADSIM  computer  program  was  developed  by  the  Con- 
vair  Aerospace  Division  of  General  Dynamics  for  Rome  Air  De- 
velopment Center  under  Contract  F-30602-73-C-0380 . It  is 
intended  to  provide  RADC  with  a simulation  which  can  be  used 
to  predict  radar  system  performance  in  many  different  envi- 
ronments. A major  feature  of  RADSIM  is  the  ability  to  in- 
clude real  world  data  for  simulation  validation  (Ref  3!xxi). 
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With  these  features,  radar  engineers  who  are  not  experienced 
simulation  designers  can  easily  model  systems  and  perform 
tradeoff  studies  without  commitment  to  expensive  hardware 
prototypes. 

The  approach  chosen  for  the  development  of  this  simula- 
tion is  to  model  the  various  components  of  a radar  system  as 
separate  modules.  The  modular  components  are  then  linked  as 
necessary  with  a driving  program  to  simulate  the  desired 
system.  Component  models  are  based  on  either  mathematical 
expressions  or  measured  data  (Ref  352-2).  RADSIM  is  coded 
in  FORTRAN  lY  and  designed  so  that  an  engineer  who  thoroughly 
understands  radar  system  analysis,  but  is  not  an  experienced 
FORTRAN  programmer  can  use  it  with  a minimum  of  assistance. 

Currently  work  is  being  done  to  extend  the  simulation 
capability  of  ECM/ECCM  performance  studies  and  to  provide 
for  interactive  operation.  The  modelling  of  Communications 
problems  has  also  been  considered.  In  1975  the  Naval  Sur- 
face Weapon  Center  (NSWC)  began  a project  to  translate  the 
program  for  use  on  a CDC  6600  computer  at  the  Dahlgren  Lab- 
oratory, but  was  forced  to  terminate  their  effort  due  to 
other  priorities.  When  RADSIM  was  selected  for  this  study, 
NSWC  provided  a partially  converted  source  program  for  use 
on  the  ASD  computer  system. 

RADSIM  Structure 

The  RADSIM  program  design  approaches  radar  modelling 
by  considering  the  radar  system  as  a collection  of  sub- 
systems divided  into  functional  elements  such  as  amplifiers. 
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filters,  and  analog  to  digital  converters  as  shown  in  Figure 
12.  Each  of  the  elements  can  be  modeled  by  specifying  its 
transfer  function  as  either  a mathematical  expression  or  as 
measured  or  derived  data  for  table  look  up.  By  linking 
these  functional  elements  to  synthesize  radar  subsystems, 
such  as  a transmitter,  receiver,  display,  or  antenna,  and 
specifying  the  appropriate  parameters  of  the  transfer  func- 
tion, a very  flexible  and  comprehensive  radar  simulation  can 
be  constructed.  In  RADSIM  each  functional  element  has  been 
implemented  as  a module  with  a distinctive  module  reference 
number  (MRN) . Usually  a module  is  a single  subroutine,  for 
instance  MRN  401  references  subroutine  NONLIN  which  repre- 
sents a zero-memory  nonlinear  device  as  shown  in  Figure  I3. 
The  transfer  function  for  this  module  is  defined  by  specify- 
ing the  corners  of  the  desired  non-linear  function  and  lin- 
ear interpolation  is  used  between  those  points. 

Data  for  signal  processing  is  contained  in  arrays 
stored  in  common  areas  for  access  by  all  modules.  All  sig- 
nal data  is  processed  in  blocks  from  these  arrays  to  reduce 
CPU  overhead  and  to  allow  snapshots  of  the  radar  signal 
states  to  be  observed  at  any  point  in  the  simulation.  Figure 
14  shows  how  several  modules  may  be  linked  to  simulate  a 
simple  radar  receiver  and  the  data  extracted  from  the  system 
at  several  points. 

There  are  two  main  programs  in  the  RADSIM  simulation. 
The  first  processes  the  input  data  which  describe  the  radar 
system  elements  and  signals  and  passes  them  to  the  simulation 
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Figure  12.  Radar  Subsystem  and  Functional  Element 


Figure  14.  Modules  Linked  for  Receiver  Model 


driver,  the  second  program.  These  programs  are  named  RADSIM 
and  KADSIM2  respectively,  but  are  stored  as  LOADR  and  SIMEX. 
(Figure  15  shows  the  flow  of  data  through  the  simulation 
activity) . 

Mathematical  Background  of  RADSIM 

Waveforms  can  be  expressed  in  real  systems  as  functions 
of  time,  S(t),  or  frequency,  S(f). 
where; 

00  -O^lTft 

S(f)=J'  S(t)e  dt  (3-1) 

_O0 

a ssimpled  waveform  of  3(t)  can  be  defined  as: 

n=oo 

S^^(t)=S(t)TE  6(t-nT)  (3-2) 

n=-oo 

where  T = the  sample  interval.  Then  the  discrete  represen- 
tation of  S(t)  over  some  time  interval  -NT<t<NT  will  be 

S^=TS(nT)  (3-3) 

The  Fourier  transform  of  this  sampled  function  S^j(t)  is 
S^j(f),  and  the  discrete  frequency  samples  will  be 

N -j2Trnl?f2t 

3^=1  TS(nT)e 
^ n=N 

-N2 

n=2 ( 3-^) 

T 

S(kQ)  is  defined  as  the  Discrete  Fourier  Transform  (DFT)  of 
S(nT)  and  H(kf))  as  the  DFT  of  h(nT)  . Discrete  convolution 
can  be  performed  on  these  two  elements  as  shown  in  Equation 
3-5.  ^ 

H(kn)S(kn)=DFT[r  S(nT)h[Mod(N, (M-n) )Tl1  (3-5) 

n=-N 
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Figure  15 . Simplified  Data  Flow  Through  Simulation 


where  MOD(N, (M-n) ) is  (M-n)  Modulo  N.  M is  the  lag  term  of 
the  convolution.  This  discrete  convolution  can  be  used  in 
determining  the  response  of  a linear  time- invariant  system  to 
some  signal  S(t)  (Ref  3S3-32).  More  details  on  the  repre- 
sentation of  specific  signal  forms  and  physical  systems  can 
be  found  in  Section  3 of  RADC  TR  76-I86.  Some  tradeoffs 
were  required  between  sampling  rates  (T)  and  signal  storage 
space  in  the  simulation.  While  higher  sampling  frequencies 
and  therefore  more  sample  points  reduces  aliasing  and  gener- 
ally leads  to  more  accurate  results,  it  was  necessary  to 
keep  the  array  sizes  within  some  reasonable  bound.  An  array 
size  of  8192  (2  ■^)  points  was  chosen.  Array  sizes  can  be 
expanded  by  allowing  overflow  into  sequential  arrays,  but 
some  simulation  flexibility  is  lost  in  the  process  (Ref  3* 

7-25). 

Ma.jor  Changes  to  RADSIM 

An  attempt  was  made  to  modify  the  program  for  the  ASD 
System  with  changes  which  would  not  significantly  alter  the 
conceptual  integrity  of  the  simulation.  It  was  necessary  to 
keep  the  changes  as  small  as  possible  and  restrict  the  scope 
of  their  effect  to  very  localized  portions  of  the  code.  This 
section  describes  the  changes. 

1)  The  HIS  FORTRAN  IV  calls  for  access  to  direct 
access  storage  have  been  changed  to  use  the  READMS/ 
WRITMS  functions  of  FORTRAN  Extended.  The  dimensioned 
areas  KEY  (6 00)  in  RADSIMl  and  RADSIM2  are  used  to 
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hold  the  table  of  direct  access  file  indexes  (Ref  9: 
1-8-9). 

2)  Honeywell  FORTRAN  IV  differs  from  GDC  FORTRAN  in 
the  use  of  subroutine  entry  point  arguments.  In 
Honeywell  FORTRAN  each  entry  point  may  have  its  own 
list  of  arguments.  In  GDC  FORTRAN  Extended,  the  entry 
point  assumes  the  arguments  of  the  parent  subroutine. 
All  subroutine  arguments  and  calling  arguments  have 
been  modified.  Where  an  argument  is  not  to  be  used  a 
dummy  variable  (DUMl  or  DUM2)  is  inserted  in  the  call. 

3)  The  Honeywell  version  of  FORTRAN  included  a func- 
tion (FLD)  for  bit  manipulation  within  words.  The  GDC 
subroutine  STRING  from  the  CDC6600  Library  was  selected 
to  perform  this  function  (Ref  10:A-14) . 

4)  A function  to  allow  the  transfer  of  data  from  a 
real  to  an  Integer  variable,  or  vice  versa,  without 
bit  position  rearrangement  was  required.  In  Honeywell 
FORTRAN  the  intrinsic  function  BOOL  performs  this  "ig- 
nore type"  requirement.  Such  a feature  is  not  avail- 
able in  GDC  FORTRAN  Extended.  It  was  provided  by 
writing  two  functions  BOOL  and  IBOOL.  They  are  doc- 
umented in  the  source  code,  and  equivalence  a location 
as  both  a real  and  integer  variable.  NSWC  changed 
some  of  the  required  function  references  to  perform 
the  function  within  the  module  code  itself.  These 
changes  are  fully  compatible  with  BOOL  and  IBOOL,  but 
could  well  be  replaced  with  function  calls  to  make  the 


program  structure  more  consistent. 

5)  File  declarations  are  made  with  control  cards  on 
the  Honeywell  system.  A PROGRAM  statement  has  been 
added  to  RADSIM  to  reflect  CDC  FORTRAN  conventions. 

6)  Honeywell  uses  a 36  hit  word.  Many  FORMAT  state- 
ments to  write  data  in  octal  or  A format  use  A6  or  06. 
Where  encountered  these  have  been  changed  to  AlO  or 
010. 

7)  The  program  uses  Namelists  to  read  input  data.  In 
some  cases  it  is  required  to  input  labels  for  graphical 
output.  Since  Hollerith  characters  are  not  allowed  in 
CDC  FORTRAN  Namelists,  the  characters  are  entered  on 
data  cards  following  the  namelist  data.  The  number  of 
words  of  Hollerith  data  in  the  label  is  read  using  the 
Namelist. 

8)  A line  printer  plotting  routine  was  added  to  aid  in 
examining  data. 

Simulation  Set-up  and  Use 

The  first  step  in  a simulation  is  to  decide  upon  a 
radar  system  or  subsystem  to  be  simulated  and  to  draw  a 
block  diagram  of  the  elements.  These  elements  must  be  spec- 
ified in  terms  of  the  functional  elements  available  in  the 
RADSIM  program.  Output  modules  are  then  added  as  desired, 
and  all  elements  are  specified  in  terms  of  their  module  re- 
ference number  (MRN).  Figure  I6  shows  the  radar  model  of 
Figure  12  in  terms  of  RADSIM  modules.  Using  the  tables  of 
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Appendicies  B and  C of  RADCTR-76-I86  the  required  parameters 


Figure  l6.  RADSIM  Modules  for  Figure  12. 


for  these  modules  can  be  generated  as  parameter  cards. 
ENDPAS,  ENDCFG,  ENDJOB  cards  are  required  at  the  end  of  the 
simulation  control  and  parameter  cards.  If  additional 
passes  or  configurations  are  run,  ENDJOB  is  used  only  after 
the  last  pass  or  configuration.  The  only  departure  from  the 
documented  instructions  for  module  parameter  cards  involves 
the  graph  labels  for  the  plot  subroutines  Modules  307-^10. 
Variables  VIL,  VDL  and  GLBL  (Independent  variable  label,  de- 
ll pendent  variable  label,  and  graph  label  respectively)  have 

; ''1 

been  replaced  with  IVIL,  IVDL,  and  IGLBL.  The  value  of 

, I these  variables  should  be  set  to  the  number  of  ten  character 

I 

. words  in  the  desired  label  (blanks  included).  The  namelist 

must  then  be  followed  by  a card  containing  a maximum  of  80 
characters  for  each  label.  IVIL,  IVDL,  and  IGLBL  and  imme- 
diately reset  to  zero  after  reading  the  label  card  so  that 
subsequent  references  to  307-310  need  not  contain  labels. 
Figure  1?  shows  an  example  of  this  usage. 

File  Structure  of  RADSIM  on  the  CPC  6600 

The  RADSIM  Program  is  stored  in  the  GDC  6600  computer 
(CSA)  under  an  ASD/EN  problem  number.  The  source  versions 
of  the  two  main  programs  are  cataloged  as  SLOADR  and  SSIMEX 
and  the  object  versions  of  these  programs  as  OLOADR  and 
OSIMEX.  Subroutines  are  stored  in  source  form  on  a file 
named  SPADLIB.  The  subroutines  may  be  located  for  modifica- 
tion in  the  Page  function  of  Editor.  There  are  two 


Libraries  of  object  subroutines.  ORADLIB,  contains  those 
subroutines  which  must  be  executed  with  OLOADR  to  produce 
the  data  file  for  execution.  ORALIBl  contains  all  the  sub- 
routines used  with  the  simulation  execution  OSIMEX.  Any  job 
could  run  using  ORALIBl/SSIMEX,  but  at  least  275K  of  core 
would  be  required.  A more  efficient  method  is  to  create  a 
new  library  containing  only  the  necessary  subroutines  for 
the  job.  The  RADSIM  documentation  suggests  some  alternate 
subroutine  groupings  which  might  offer  advantages. 

Core  and  Time  Requirements 

An  average  RADSIM  job  requires  approximately  225000g 
words  of  central  memory  and  35-40  seconds  per  pass  through 
the  simulation.  To  permit  more  accurate  estimates,  calls  to 
SECOND,  and  TIME  are  included  in  the  SIMEX  program.  Use  of 
a partial  load  map  can  also  assist  in  estimating  additional 
memory  requirements  for  job  extensions.  System  usage  can  be 
reduced  by  executing  data  decks  with  the  loader  (OLOADR  and 
ORADLIB)  only  to  verify  their  correctness.  In  this  way  er- 
rors in  variable  names  in  the  namelists  can  be  found  and  cor- 
rected prior  to  full  scale  execution.  About  70K  of  core  and 
less  than  ten  seconds  of  CPU  time  should  be  sufficient  for 
normal  jobs  of  this  nature.  Experience  with  the  ASD  com- 
puter center  operation  indicates  that  normally  a full  RADSIM 
job  will  only  execute  overnight  unless  special  priorities 
are  obtained. 


Suggested  Changes  to  RADSIM 


The  following  are  changes  which  have  not  "been  per- 
formed hut  would  inhance  simulation  integrity  or  operation, 
a)  Convert  all  in-MODULE  code  "ignore- type"  code  to 
calls  to  BOOL  or  IBOOL. 

h)  Check  out  all  remaining  FOFUVIATS  for  incorrect  word 
size  in  0 or  A FORMATS. 

c)  Consider  some  form  of  segmentation  or  overlaying 
to  reduce  core  requirements.  This  may  he  complicat- 
ed hy  large  common  area  requirements. 

d)  Consider  taking  advantage  of  the  increased  pre- 
cision of  the  60  hit  word  over  the  36  hit  word. 

Verification  of  Correct  Performance 

Since  the  purpose  of  this  thesis  included  the  simula- 
tion of  synchronized  pseudo-random  noise  against  a DICKE-FIX 
radar  receiver,  the  DICKE-FIX  Sample  Problem  of  Volume  I 
Part  3 was  used  as  verification  of  the  simulations  correct 
performance.  Figures  I8  and  19  show  the  comparative  results 
on  the  Honeywell  and  CDC  computers.  The  input  values  were 
the  same,  and  the  outputs  were  very  close.  The  data  from 
the  CDC  was  hand  plotted  while  the  HIS  data  was  plotted  on 
an  HP9820  calculator  and  x-y  plotter.  Since  this  is  essen- 
tially what  will  he  run  for  the  experiment,  the  simulation 
was  deemed  correct  for  this  case. 


Figure  18.  RADSIM  Response 


Figure  19 • Response  of  RADSIM  to  TGT  Pulse 


IV.  Experiment  Design 


The  second  major  task  of  this  study  is  to  examine  the 
effectiveness  of  synchronized  pseudo-random  noise  jamming 
when  employed  against  a radar  equipped  with  a simple  ECCM 
device,  the  DICKE-FIX.  To  do  this  an  experiment  will  he 
performed  using  the  RADSIM  program.  This  section  will  de- 
scribe the  design  of  the  experiment.  The  five  phases  to  be 
covered  are: 

1)  Selection  of  the  radar  parameters. 

2)  Selection  of  radar  model  parameters. 

3)  Design  of  a synchronized  pseudo-random  noise  jammer. 

4)  Set-up  of  the  RADSIM  model. 

5)  Specification  of  the  experimental  procedure. 


Radar  Parameters 

The  parameters  of  the  victim  radar  for  this  experiment 
are  shown  in  Table  II.  Only  the  receiver  and  display  sub- 
systems are  to  be  modeled.  Since  the  RF  waveform  carries  no 
significant  information  it  will  be  neglected  (Ref  2:208). 


Table  II 

Radar  Parameters 
Intermediate  Frequency  (IF) 

IF  Bandwidth 

Pulse  Repetition  Frequency 
Pulse  Width 

Integration  Efficiency 


20  MHZ 
1 MHZ 
300  PPS 

1200  Nano-seconds 
0.95 
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The  intermediate  frequency  and  bandwidth  are  chosen 
to  be  representative  of  search  radars  which  might  be  en- 
countered in  USAF  operations.  They  do  not  correspond  to 
any  specific  system.  The  pulse  width  was  set  using  the  re- 
lationship (Ref  7i47)« 

T = 1.2/B  = 1200  ns.  (4-1) 

where:  B = radar  bandwidth 

T = radar  pulse  width 

With  the  basic  "radar  parameters  of  interest  specified, 
the  receiver  system  block  diagram  can  be  drawn.  Using  the 
model  of  Table  II,  the  radar  receiver  subsystem  is  shown  in 
Figure  20.  The  target  return  signal,  thermal  random  noise 
and  jamming  are  summed  at  the  input  of  the  receiver.  The 
signal  then  passes  through  the  wideband  limiter  (bandwidth  > 
10  MHZ),  is  filtered  at  the  intermediate  frequency  to  recover 
the  target  signal,  and  then  is  detected  using  a full  wave  de- 
tector. The  energy  in  the  signal  is  integrated  to  form  a 
single  pulse  waveform  and  amplified  by  some  arbitrary  gain. 
This  amplified  signal  is  then  applied  to  a post-detection 
video  sweep  integrator  which  integrates  eight  pulses  for  the 
display.  Four  cases  will  be  considered  in  this  experiment: 

1)  Thermal  random  noise  only 

2)  Target  plus  thermal  random  noise 

3)  Thermal  random  noise  plus  synchronized  pseudo- 
random noise  jamming 

4)  Target  plus  thermal  random  noise  plus  synchronized 


Figure  20.  Radar  Receiver  Model 
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pseudo-random  noise  jamming 
Radar  Model  Parameters 

To  implement  the  radar  model,  the  parameters  of  each 
element  must  he  specified.  Some  of  the  elements,  the  wide- 
band limiter,  the  target  signal  shape,  and  the  IF  filter 
parameters  are  characteristic  of  the  radar  and  will  remain 
unchanged  in  any  iteration  of  the  experiment.  Others  such 
as  the  linear  amplifier  gain  and  the  integration  efficien- 
cies of  the  two  integrators  may  be  varied  but  are  not  crit- 
ical parameters  for  this  investigation  and  will  remain  fixed 
throughout  the  experiment.  The  peak  voltage  of  the  target 
signal  and  the  mean  and  variance  of  the  noise,  both  thermal 
and  SPRN  are  parameters  which  may  be  adjusted  to  provide  the 
desired  signal  to  noise  ratios  (SNR)  for  parametric  analysis. 
Since  the  synchronized  pseudo-random  noise  will  be  generated 
off-line  to  the  simulation  its  characteristics  will  remain 
fixed. 

The  Wideband  Limiter.  The  wideband  limiter  character- 
istic is  that  of  a hard  limiter.  All  positive  input  voltages 
become  and  all  negative  voltages  -V^  (Ref  7tl20).  Figure 
21  shows  the  voltage  transfer  function  of  the  limiter.  Its 
frequency  characteristic  is  assumed  linear  over  the  band- 
width of  the  simulation. 

The  Intermediate  Frequency  (IF)  Filter.  The  IF  filter 
was  selected  to  be  a fourth-order  Butterworth  filter.  This 
form  is  maximally  flat  in  the  passband  and  closely  approxi- 
mates the  envelope  of  the  power  spectral  density  of  the 


rectangular  pulse  transmitted  by  the  radar  (Ref  11;  71). 
This  can  be  shown  by  taking  the  Fourier  Transform  of  the 


1200  ns  rectangular  unit  pulse. 

S(t)=p[t,  t-1200] 
S(f)=1200  sine  (1200f) 


(^-2) 

(^-3) 


S(f),  shown  in  Figure  22,  closely  approximates  the  form  of 

the  4th  order  Butterworth  filter  shown  in  Figure  23-  (Ref 

11:71).  The  lowpass  prototype  of  this  filter  is  specified 

by  a transfer  function  H(s). 

H(S)=  1 _ 

sVaS^tS^+cS+l 

where  a,  b,  and  c are  coefficients  which  may  be  calculated 
or  found  in  tables.  In  this  case,  the  values  from  tables  of 
Reference  11  were  used.  The  prototype  low  pass  filter  is 
transformed  to  a bandpass  filter  by  the  frequency  trans- 


formation! 


S=s®+(i)(,® 


(4-5) 


S=J(i3 

a)o=the  center  frequency  of  the  bandpass  filter 
B=the  bandwidth  of  the  filter  in  radian  frequency 
The  transformed  filter  response  is  then 


( s^+Wo ^ ) +a(  s^+a)o ^ ) ^+b(  s^+Wo ^ ) ^+C ( s^+cOo ^ )+l 


s +As  +Bs  +Cs^+Ds  +Es'^+Fs  +u)e 
For  simulation  purposes  these  coefficients  are  cal- 
culated in  terms  of  frequency  rather  than  radian  frequency 


1 

Q 

and  the  units  are  in  Gigahertz  (10^).  Table  III  lists  the 
coefficients . 


% Table  III 

^ IF  Filter  Transfer  Function  Coefficients 

A=ab=0 . 26131259X10"^ 

B=bB^+4fo^=0. 16034142x10"^ 

C=CB^+3afo^B=0. 31383642x10"-^ 

I D=2B^fo^+6fo^+B^=0. 96273237x10"^ 

E=CB^fo^+3aBf 0^=0. 12553457x10"® 
F=bB^fo^+4fo^+0.25654627X10"9 
j G=abfo^=0. 16724006x10"^^ 

I fo®=0.256000X10"^3 

B^=0. 10000000X10"^^ 


Factoring  Equation  (2-6b)  with  the  AFIT  PROGRAM  POLY,  pro- 
vided the  complex  poles  for  the  IF  filter.  Using  another 
AFIT  PROGRAM,  FREQR,  the  filter  was  verified  to  have  es- 
sentially unity  gain  at  20  MHZ  with  _^3db  points  at  19.5MHZ 
and  2O.5MHZ  respectively.  Figure  24  shows  the  filter  trans- 
fer function  plotted  in  the  range  15-25MHZ . 

Integration  Efficiency.  Integration  efficiency  will 
be  set  at  0.95* 

Target  and  Noise  Signal . The  target  signal  is  speci- 
fied to  be  a rectangular  pulse  at  the  20  MHZ  center  fre- 
quency of  1200  Nanoseconds  duration.  The  rise  time  is  speci- 


(ZHM) 


the  IF  filter  will  provide  the  desired  signal  to  noise 
ratio. 


The  Synchronized  Pseudo-Random  Noise  Jammer 

By  the  theory  developed  in  Section  II,  there  are  two 
major  constraints  which  must  be  met  effective  pseudo-random 
noise  jamming.  A bandwidth  constraint  which  defines  the 
point  of  the  first  null  of  the  jamming  noise-power  spectral 
density  function  and  a synchronization  constraint  that  is 
imposed  by  the  requirement  that  the  period  of  the  sequence 
equal  the  pulse  repitition  interval  of  the  radar.  These  can 
be  expressed  as  a pair  of  simultaneous  equations  in  terms  of 
the  clock  frequency  of  the  shift  register  (Ref  2:138). 

(^-7) 


f =kBr 
2 


f =2^f 


f =clock  frequency 


B=radar  bandwidth 

r=number  of  stages  in  the  shift  register 
f^=pulse  repitition  frequency  of  the  radar 
k-a  constant  to  increase  the  bandwidth  of  the  jammer 
by  some  factor  to  allow  for  inaccuracies  in  deter- 
mining the  radar  parameters  (set-on  constant) 
Figure  25  shows  the  results  of  plotting  the  values  of  these 
two  equations  for  the  radar  parameters  of  Table  II.  Assum- 
ing that  k=10  is  a reasonable  value  for  the  set-on  constant, 
the  clock  frequency  and  shift  register  length  are  ''92MHZ  and 
18  respectively.  By  choosing  the  primitive  irreducible 
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j Figure  25.  Plot  of  Synchronization  And  Bandwidth 

* r!nnc!+.T’a  1 n+  'P'niisi +.  i nn  1 1 1 -h i nn<3 


polynomial  of  degree  18  with  the  smallest  number  of  non-zero 
terms  from  the  Tables  of  ENAD  TR-7^-01  (Ref  5 ••97)  the  char- 
acteristic function  is  found  to  bej 

f(x)=X^®+X^^+l  (4-9) 

This  leads  to  a shift  register  jammer  of  the  form  shown  in 
Figure  6.  This  jammer  was  implemented  in  a computer  pro- 
gram to  provide  8192  samples  values  of  the  jamming  signal  in 
a 4000  observation  interval  around  the  received  pulse. 
Appendix  A describes  to  computer  program  PRN  which  calcu- 
lates the  sample  values  and  stores  them  in  a form  suitable 
for  RADSIM. 

The  RADSIM  Model 

Figure  26  shows  the  radar  model  of  Table  II  specified 
in  terms  of  RADSIM  Module  Reference  Numbers.  Table  IV  de- 
scribes the  modules  and  their  specific  functions.  It  also 
specifies  the  sequence  in  which  the  modules  are  to  be  exe- 
cuted. Appropriate  output  modules  have  been  added  to  allow 
examination  of  the  wave  forms  at  two  points  in  the  simula- 
tion. Table  V lists  the  control  cards  and  simulation  cards 
necessary  to  perform  the  simulation.  Figure  27  shows  the 
data  flow  for  simulation  execution.  Note  that  the  sampled 
data  for  the  jamming  wave  form  has  been  stored  on  a tempo- 
rary file  (TAPE  3) I and  is  rewound  and  repeated  every  pulse 
repetition  interval. 

Experimental  Procedure 


The  experiment  will  consist  of  two  phases.  The 
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Figure  26.  Radar  Model  as  run 


t 


‘ i 

I 


I ' 


Table  IV 


Module  Specifications  and  Sequence 

Seq.  No. 

Module  Name/MRN 

Purpose/Specification 

1. 

system/301 

Initialize  System  Parameters 

2. 

RRAND/102 

Initialize  Random  Number  Gener- 

ator 

3. 

FGENXY/402 

Target  Signal  Generator 

4. 

ADDRND/235 

Add  Random  Noise  to  Target 

Signal 

5. 

REinflND  3 

Rewind  PRN  Jammer  Sample  File 

6. 

STOREA  3 

Store  Jammer  Samples  in  Array 

XA 

7. 

ADDA/207 

Add  TGT  + Noise  + Jammer 

8. 

NONLIN/401 

Wideband  Limiter 

9. 

ERGYCP/116 

Measure  Input  Energy 

10. 

ZFFT/202 

Transform  to  Frequency  Domain 

11. 

FILT/407 

4th  Order  Butterworth  Filter 

IMHZ  BW  @ ZOMHZ 

12. 

ZIFFT/203 

Transform  to  Time  Domain 

13. 

ERGYCP/116 

Measure  Output  Energy 

14. 

FWDET/416 

Full  Wave  Envelope  Detector 

15. 

INGTOR/409 

Integrator  (-95  Efficiency) 

16. 

LAMPCD/456 

Linear  Amplifier 

17. 

* PLOTTR/307 

Plot  Presv/eep  Integration 

Signal 

18. 

* CUMDIS/208 

Generate  Histogram  of  Signal 

Voltages 

J 
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Seq.  No. 

Module  Name/MRN 

Purpose/Specification 

19. 

* OUTCUM/211 

Generate  Cumulative  Distribu- 
tion of  Signal 

20. 

* PLOTTR/309 

Plot  CDF 

21. 

ATOD/216 

Analog  To  Digital  Conversion 

22. 

DTOA/218 

D to  A (Set  Saturation  Level, 

Reduce  Number  of  Points) 

23. 

NCSWPl/436 

Post  Detection  Video  Sweep  In- 
tegrator 

24. 

* PLOTTR/309 

Plot  Integrated  Signal 

25. 

* CUMDIS/208 

Generate  Histogram  of  Inte- 
grated Signal 

26. 

* OUTCUM/211 

Generate  CDF  of  Integrated 

Signal 

27. 

* PLOTTR/309 

Plot  CDF 

* 

Omitted  on  Passes 

1-7  to  reduce  output  requirements 

! 
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Table  V 

Control  Cards  for  RADSIM  Execution 


JOB  CARD 

FTN.  Used  to  Execute  SPRN  Generator 

LGO.  Not  used  W/WGN  Cases 

ATTACH  OLOADR 
ATTACH  ORADLIB 

ATTACH  CC6600,  ID=X654321,  SN=ASD 

LIBRARY  ORADLIB,  CC66OO 

OSIMEX 

7/8/9/ 

SPRN  Generator  Source  And  Assembly  Language  Deck 
7/8/9 

DATA  Card  for  SPRN  Generator 
7/8/9 

Simulation  Control  Card  Deck 


S/7/8/9 


Figure  27-  Data  Flow  as  run 
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first  phase  will  be  calibration.  The  output  of  the  IF  fil- 
ter (MRN203)  will  be  measured  for  two  cases,  White  Gaussian 
Noise  (WGN)  and  WGN  + Synchronized  Pseudo-Random  Noise.  By 
calculating  the  energy  density  for  each  of  these  cases,  the 
desired  signal  to  noise  ratio  may  be  obtained  by  adjusting 
the  target  return  signal  voltage.  A signal  to  noise  ratio 
of  2 will  be  studied  since  it  should  show  the  largest  change 
in  performance  (Ref  2t273)*  The  final  phase  will  be  experi- 
mentation in  which  the  simulation  will  be  exercised  for  the 
cases  shown  on  Table  VI. 


Table  VI 

Simulation  Cases 


1. 

WGN 

(No  DICKE-FIX) 

2. 

WGN 

(With  DICKE-FIX) 

3. 

WGN 

+ SPRN  (No  DICKE-FIX) 

4. 

WGN 

+ SPRN  ( With  DICKE-FIX) 

5. 

WGN 

+ TGT  (No  DICKE-FIX)  (SNR=2) 

6. 

WGN 

+ TGT  (With  DICKE-FIX)  (SNR=2) 

7. 

WGN 

+ TGT  + SPRN  (No  DICKE-FIX)  (SNR=2) 

8. 

WGN 

+ TGT  + SPRN  (With  DICKE-FIX)  (SNR=2) 

Using  the  data  gathered  from  this  sequence  of  simula- 


tions, a set  of  receiver  operating  characteristic  (ROC) 
curves  will  be  constructed  for  cases  five  through  eight. 
These  ROC's  will  be  used  to  verify  the  theory  and  demon- 
strate the  effectiveness  of  synchronized  pseudo-random  noise 
jamming  against  the  DICKE-FIX  receiver. 
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V.  The  Experiment,  Data  Analysis , and  Results 


In  this  section  the  conduct  of  the  experiment  will 
be  explained  and  the  data  obtained  analyzed.  There  were 
two  major  phases  to  the  experiment,  as  described  in  Section 
IV.  The  first  consisted  of  a series  of  measurements  to  de- 
termine the  amount  of  energy  present  in  the  White  Gaussian 
Noise  (WGN)  and  Synchronized  Pseudo-Random  Noise  (SPRN)  sig- 
nals so  that  an  appropriate  target  signal  amplitude  could  be 
established,  and  related  to  the  signal  to  noise  (SNR)  ratio 
of  Dr.  Berrie's  experiments.  The  second  phase  consisted  of 
running  each  of  the  cases  described  on  “liable  VI  of  Section 
IV  with  the  appropriate  target  signal  amplitude.  Simula- 
tion output  data  consisted  of  the  cumulative  probability 
distribution  of  the  pre  and  post -integration  signal  voltages 
for  the  eighth  pulse.  The  data  was  recorded  in  both  graph- 
ical and  point  list  form.  This  data  was  converted  to  P^^^^ 
and  P^  curves  similar  to  Figure  11,  and  an  ROC  was  drawn  for 
each  case.  These  ROC's  were  then  compared  to  the  data  ob- 
tained in  Berrie ' s experiments  and  related  to  the  analyt- 
ical values  he  established. 

Determination  of  WGN  Energy  Levels 

Reference  3 provides  a formula  for  establishing  the 
energy  of  WGN  noise  generated  by  the  simulation. 
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(5-1) 


E = (2^^  ) (0^)  (1/FS)  (Watt-Nanoseconds) 
where:  E = total  energy 

2 = no.  of  samples  as  determined  by  input  param- 

eters (2'^^=2^^=8192 
0®  = noise  variance 
FS  = sample  rate 

Energy  density  can  then  be  determined  using  the  energy  ob- 
tained from  Equation  (5-1)* 

N = E/FS  ( Watt-Nanosec onds/GHZ)  (5-2) 

N = energy  density 

The  value  selected  for  the  noise  variance  was  0.5,  the  sam- 
ple rate  was  2.0  GHZ,  and  N®  was  13*  Thus  by  (5-1)  and 
(5-2),  including  the  factor  (2400/4096)  to  account  for  the 
sampling  of  energy  during  the  target  pulse  only; 

E = (2^^) (0.5*) (1/2) (2400/4096)  = 600  W-NS  (5-3a) 

N = 600/2  = 300  W-NS/GHZ  = 0.3W-NS  (5-313) 

MHZ 

Measurement  of  ten  independent  runs  of  WGN,  taking  the 
energy  of  the  signal  at  the  input  and  output  of  the  IF 
filter,  provided  the  data  of  Table  VII.  Each  run  consisted 
of  eight  passes  through  the  simulation  with  only  the  mean 
and  standard  deviation  of  the  energy's  obtained  shown. 


Table  VII 

Statistics  of  Ten  Independent  White  Gaussian  Noise 


Sample  Sets 

Into  and  Out  ( 

of  the  IF 

Filter 

IDENT 

^IN 

(Ei)W-NS 

7 

Eout(^o)W-NS 

MEAN 

STD  DEV 

MEAN 

STD  DEV 

FMA 

593.^0 

25.17 

0.4l 

0.24 

FMB 

602.8? 

15.22 

0.70 

0.34 

EMC 

599.03 

9.92 

0.65 

0.51 

FMD 

588.08 

13.79 

0.67 

0.31 

FME 

595.54 

8.95 

0.37 

0.28 

FMF 

602.73 

16.94 

0.62 

0.23 

FMG 

600.57 

11.35 

0.71 

0.27 

FMH 

602.91 

12.53 

0.69 

0.44 

FMI 

601.85 

12.73 

0.49 

0.20 

FMJ 

594.46 

13.47 

0.25 

0.52 

Calculating  the  mean  and 

standard 

deviation  of 

these 

samples  gave: 

- 598.14 

s = 5.06 

(5-4a) 

Eo  = 0.56 

s = 0.16 

(5-4b) 

n = 0.29907 

Watt-Nanoseconds/MHZ 

(5-4c) 

These  values  of  and  n compared  favorably  with  the  theo- 
retical value  of  Equations  (5-3a.  5-3b) . The  noise  band- 
width (Bj^)  of  the  receiver  model  was  then  computed  using 
Equation  (5-5)* 


6^ 


i 


= Eo/n  {5-5) 

= 0.56/0.30 

= 1.87MHZ 

Since  the  signal  energy  and  noise  energy  are  measured 
over  the  same  period,  the  ratio  of  signal  power  to  noise 
power  (SNR)  is  equivalent  to  the  ratio  of  signal  energy  to 
noise  energy.  For  a noise  bandwidth  of  I.87  MHZ,  the  energy 
ratio  at  the  input  to  the  IF  filter,  equivalent  to  a post 


IF  filter  SNR  of  2,  is 


= Z 

s s 


Eo  nB,, 


Eg  = (2)(1.87)(0.3)  = (2)(1.87)(0.3) 
= 1.122  Watt  -Nanoseconds 


(5-6a) 


(5-6b) 


Signal  amplitude  was  then  calculated  using  Equation 


(5-7) (Ref  14:321). 


Eg  = 1/2  AT 


(5-7) 


A = signal  amplitude 
T = signal  deviation 

For  E = 1.122  Watt  -Nanoseconds  and  T = 1200  Nanoseconds 
s 

(from  Table  II )i 


A = V2 {E, 


(5-8) 


= /"2  (1.122)  = 0.043  volts 

1200 

The  value  of  A from  Equation  (5-8)  was  then  used  as  VPEAK 
for  a SNR  = 2 in  White  Gaussian  Noise. 

Since  the  amount  of  energy  within  the  IF  filter  band- 


T 


width  varied  from  run  to  run,  a high  value,  low  value  and 
near  mean  value  were  chosen  to  be  run  for  the  experiment. 

The  high  was  obtained  from  run  FMG,  the  low  from  FMJ  and 
the  closest  to  the  mean  from  FMF. 

Determination  of  Synchronized  Pseudo-Random  Noise  Energy 
Levels 

To  determine  the  energy  present  in  the  SPRN  jamming 
signal  a calculation  was  performed  using  a formula  from  the 
users  manual  for  the  Hewlett  Packard  3722A  Noise  Generator 
(Ref  1353-7) • Its  characteristics  are  very  similar  to  the 
noise  generator  used  in  this  experiment.  The  Noise  Power 
Density  of  the  Pseudo-Random  Noise  generator  is  given  by 
Equation  (5-9). 

NPD  = 2a®  T (5-9a) 

NPD  = noise  power  density  in  V^/HZ 
a = the  shift  register  amplitude 
T = the  clocking  rate  of  the  shift  register 

Thus : 

NPD  = (2)(1)(V92X10®)  = 2.174  X 10“^  V®/HZ  (5-9b) 
Energy  density  is  the  product  of  the  noise  power  density 
times  the  interval  over  which  the  measurement  is  taken,  so 
by  Equation  (5-10) 

E = (NPD) (Time)  (5-10) 

= (2.174  X 10"^  Watts) (2400  N.S.) 

= 521.8  Watt  - Nanoseconds 

To  measure  the  actual  output,  the  noise  generator  was 
cycled  through  a full  period  and  the  register  state  sampled 
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every  37?  shifts.  Ten  of  these  states  were  selected  at 
random  and  used  as  the  starting  point  for  a simulation. 

The  average  energy  for  each  of  these  runs  and  the  standard 
deviation  is  shown  in  Table  VIII. 


Table  VIII 

Statistics  of  Ten  Independent 
Pseudo-Random  Noise  Jamming  Signals 


IDENT 

MEAN 

Ejj^(Ei)W-NS 

STD  DEV 

F 

^OUT 

MEAN 

(Eo)W-NS 

STD  DEV 

FMK 

771.64 

18.3^ 

25.74 

2.45 

707.5^ 

14.30 

8.36 

2.10 

FMM 

689-50 

9.64 

7.95 

1.07 

679.90 

12.70 

7 . 66 

1.78 

691.01 

14.63 

7.03 

1.40 

FMP 

692.76 

11.13 

10.79 

1.90 

FMQ 

731.72 

13.28 

10.51 

1.52 

FMR 

669 . 68 

10.45 

5.28 

0.94 

FMS 

690.29 

12.86 

55.08 

4.28 

FMT 

691.98 

14. 84  . 

11.21 

2.30 

Averaging  these  samples  provided  a mean  energy  into 
the  IF  filter  of  701.59  Watt-  Nanoseconds  and  a mean  energy 
out  of  14-.96  Watt-  Nanoseconds.  Standard  deviations  were 
29.62  Watts-  Nanoseconds  and  15.20  Watt  - Nanoseconds  re- 
spectively. These  values  of  means  did  not  compare  closely 
to  the  calculated  value  of  energy  in  the  SPRN  signal  from 
Equation  (5-9).  Equation  (5-9)  however  assumed  that  the  en- 
tire sequence  was  generated.  Berrie  showed  that  for  short 
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subsequences  of  maximal  length  sequences,  the  Chi-square 
value  can  be  very  high  (Ref  2tl09)»  One  would  therefore, 
not  expect  any  individual  sequence  to  demonstrate  an  exact 
correspondence  to  the  full  sequence . Further  verification 
of  this  property  was  shown  by  taking  the  average  value  of 
the  modulating  signal  for  one  of  the  runs.  For  a true  max- 
imal length  sequence  the  average  value  should  be  very  close 
to  zero.  Examination  of  the  average  values  of  the  subse- 
quences showed  a distinct  bias  which  would  add  some  energy 
to  the  waveform  that  would  not  exist  for  a true  maximal 
length  sequence.  There  were  695  possible  independent  subse- 
quences of  which  only  ten  were  chosen.  It  is  quite  possible 
that  examination  of  more  subsequences  would  produce  an  aver- 
age value  closer  to  the  computed  value. 

The  target  signal  amplitude  was  set  in  the  same  manner 

as  for  White  Gaussian  Noise,  using  the  calculated  values  for 

energy  and  energy  density.  The  bandwidth  (B.)  of  the  jammer 

(J 

was  10.22  MHZ.  Thus: 

’^SPN  ^ ®/®j  (5-11) 

= 521.8/10.22 

= 51.1  Watt  -NS 
MHZ 


Substituting  into  Equation  (5-6a) 

' "SPN^N^ 

= (2)(1.87)(51.1) 

= 191.15  Watt  - NS 

And  computing  the  amplitude  A as  in  Equation  (5-8) 


(5-12) 
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^ 

A = /I  {191.15  = 0.56  volts  (5-13) 

1200 

VPEAK  for  the  Synchronized  Pseudo -Random  Noise  cases  was  set 
at  0.56  volts. 

Table  VIII  showed  that  the  amount  of  energy  in  the  IF 
filter  bandwidth  varied  greatly  and  no  mean  value  was  ever 
generated.  Two  samples  were  used.  FMK  and  FMT  bracketed 
the  calculated  mean  value . 


Conduct  of  the  Experiment 

Each  case  was  run  using  the  parameters  which  defined 
the  energy  levels  and  the  appropriate  target  signal  voltage. 
Limits  on  the  output  data  histogram  elements  were  changed  as 
necessary.  Each  case  from  Table  VI  was  generated. 

Results  and  Analysis  of  the  Data 

The  data  from  the  simulation  consisted  of  a point 
list  and  graph  describing  the  cumulative  probability  distri- 
bution (c.p.d.)  of  the  voltage  levels  both  into  and  out  of 
the  pulse  integrator  for  the  eighth  pulse  integration.  The 
c.p.d.  indicated  the  probability  that  any  voltage  would  be 
less  than  a specified  voltage,  as  shown  in  Equation  (5-1^) 
and  Figure  9. 


P(v<V)  *=  / f(v)  dv 


(5-1^) 


Probability  of  detection  or  false  alarm  as  defined  in  Equa- 
tions (2-23)  and  (2-24)  may  be  calculated  from  the  Simula- 
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tion  c.p.d.  by  simply  subtracting  any  value  on  the  c.p.d. 


from  1.  A pair  of  curves  representing  and  were  de- 


veloped for  each  sample  case  as  shown  in  Figure  11.  These 


curves  were  then  used  to  generate  a Receiver  Operating  Char- 


acteristic Curve  that  combined  the  pre  and  post -integration 


performance  of  the  individual  sample  case.  For  WGN  there 


was  a high,  low,  and  near  mean  sample  for  each  case.  In  the 


SPRN  cases  there  was  a high  mean  and  a low  mean  sample. 


These  samples  were  averaged  to  produce  an  average  Receiver 


Operating  Characteristic  curve  for  each  of  the  cases  five 


through  nine  of  Table  VI.  These  ROC’s  are  shown  in  Figures 


28  through  31*  Each  performance  curve  was  also  compared  to 


Berrie's  data  to  evaluate  the  correspondence  between  the 


RADSIM  simulation  and  the  values  obtained  from  theory  and 


analog  experiments.  Figures  32  through  35  show  these  com- 


parisons. Performance  comparisons  for  the  cases  with  and 


without  the  wideband  limiter  are  made  in  Figures  3^  through 


Discussion  of  Results 


A true  ROC  should  go  to  zero  in  its  lower  left  hand 


corner  indicating  P^^^  “ ^d  ~ non-zero  value  of  the 


observed  ROC's  is  due  to  truncation  of  simulation  values  to 


fit  into  discrete  histogram  elements,  of  which  there  are  only 


a finite  number.  Better  data  could  be  obtained  in  these 


ranges  by  selecting  different  histogram  parameters.  All  of 


RADSIM  ROC's  exhibited  this  characteristic. 


Figures  32  through  35  show  that  the  RADSIM  results 
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Probability  of  Detection 


Figure  28.  Average  Periorrr.ance  in  WGN 
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Probability  of  Detection  (P 


Probability  of  False  Alarm  (P^g^) 


y I-  1 < 2'>.  Average  Performance  in  Synchronized 

Pccudo-Pandom  Noise 
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Probability  of  Detection  (P 


0.8  I 


0.6  I 


Probability  of  False  Alarm  (P^g^) 


Probability  of  Detection  (P 


Probability  of  False  Alarm  (P^^) 


\ Figure  32.  Comparison  of  Average  Post-Integration 

5 Performance,  in  WGN.  RADSIM  vs. 

ANALOG  Simulation. 


Probability  of  Detection 
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Probabili 


Figure  Comparison  of  Average  Post-Integration 

Performance  in  SPRN.  RADSIM  vs.  ANALOG 
Simulation 
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II 


Pro'bability  of  Detection  (P 


Probability  of  False  Alarm  (P^g^) 


Figure  35*  Comparison  of  the  Average  Pre-Integration 
j Performance  in  SPRN.  RADSIM  vs.  ANALOG 

■i  Simulation 


Figure  3^.  Comparison  of  Average  Post-Integration 
Performance  in  WON  VJith  and  Without 
Limiting.  SNR  = 2 
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Figure  3?.  Comparison  of  Post-Integration  Performance, 
I With  and  VJithout  Limiting,  in  Synchronized 
H Pseudo-Random  Noise. 


i'ii 


'5! 


■(  J 


were  generally  compatible  with  the  theoretical  development 
of  Section  II.  An  improvement  in  performance  is  evident 
when  the  receiver  is  operating  in  WON,  but  performance  de- 
grades significantly  in  SPRN  for  approximately  the  same  SNR. 

The  comparisons  of  performance  between  RADSIM  results 
and  the  theoretical  data  developed  by  Berrie  show  a good 
correspondence  for  WGN  and  the  same  tendencies  as  Berrie  for 
SPRN.  The  actual  numerical  differences  are  probably  the  re- 
sult of  differences  in  the  actual  signal  to  noise  ratios 
used . 

The  use  of  a limiter  did  not  demonstrate  any  perform- 
ance improvement  over  the  non  limiting  receiver.  In  the  WGN 
cases  a moderate  performance  degradation  occurred  as  implied 
by  Nathanson  (Ref  7:122).  In  the  SPRN  case  a very  slight 
degradation  was  noted. 
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VI.  Conclusions  and  Recommendations 


The  results  described  in  Section  V show  that  the 
original  objectives  of  this  thesis  have  been  met.  The 
RADSIM  program  has  been  successfully  adapted  to  the  ASD  CDC 
6600  computer  system  and  shown  to  be  a useful  tool  for  the 
analysis  of  a search  radar  receiver.  The  results  of  Doctor 
Berrie's  investigations  of  the  effect  of  Pseudo-Random  Noise 
on  radar  performance  have  been  approximately  duplicated,  al- 
though an  exact  numerical  correspondence  was  neither  ex- 
pected nor  achieved.  The  use  of  a DICKE-FIX  ECCM  device  has 
been  shown  to  be  ineffective  against  synchronized  Pseudo- 
Random  Noise  Jamming  of  the  form  employed  in  this  experiment. 

Further  studies  are  required  to  demonstrate  the  capa- 
bilities of  SPRN  jamming  against  other  ECCM  devices.  Also, 
the  spectral  characteristics  of  subsequences  of  maximal 
length  sequences  should  be  investigated.  The  amount  of  en- 
ergy in  the  bandwidth  of  the  IF  filter  fluctuated  greatly. 

In  a range-gated  radar  receiver  the  target  would  appear  in 
some  range  bins  and  be  masked  in  others . Criteria  should  be 
developed  for  providing  the  masking  in  the  range  bins  con- 
taining the  target. 
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Appendix  A 


This  appendix  describes  the  program  SPRNGEN  used  to 
provide  samples  of  pseudo-random  noise  jamming  for  use  in 
the  RADSIM  program. 

Program  Structure 

SPRNGEN  consists  of  a main  program  coded  in  FORTRAN, 
and  an  assembly  language  subroutine.  The  main  program  shown 
in  Figure  38,  reads  the  data  which  provides  the  parameters 
for  the  generator,  sets  up  the  required  constants  and  calls 
the  assembly  language  subroutine  whenever  a register  shift 
is  required.  The  assembly  language  subroutine,  shown  in 
Figure  39  shifts  the  register  one  bit  and  returns  with  the 
register  sum  (NCNT)  and  a bit  of  the  pseudo-random  code 
(NBIT).  This  sum  is  used  to  provide  DSB  modulation  to  a 
sampled  carrier  signal.  Samples  are  written  to  TAPE3f  and 
a count  of  the  sums  of  the  register  states  maintained. 

Input  Data 

Data  is  read  by  a single  list  directed  READ  statement. 
Table  IX  lists  the  input  data  variables  and  their  functions 
and  restrictions,  if  any. 
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Table  IX 

SPRNGEN  Input  Variables 

CLKRT  - The  clock  rate  of  the  shift  register. 

SAMPRQ  - The  sampling  frequency  of  the  output  signal.  It 

must  be  equal  to  the  simulation  sample  frequency  FS. 

CARPRQ  - The  carrier  frequency  of  the  jamming  signal. 

N2 

NSAMP  - The  number  of  samples  desired.  It  must  equal  2 
where  N2  is  the  same  as  the  simulation  variable. 

N - The  number  of  stages  in  the  shift  register. 

MASK  - The  Feedback  Mask  of  the  shift  register.  It  should 

be  the  decimal  equivalent  of  a binary  one  in  the 
position  of  each  feedback  tap. 

NSEQ  - The  initial  shift  register  state  as  a decimal 
number.  It  may  not  be  zero. 

TORIG  - Not  used. 
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IBIS  PAQS  IS  BSST  QUALITY  FRACtlC^W 
F»ailCOPYFURN.ISHBI)XOQD.O  - 


PS'OGF  ftM  S^RNGE  M(  IN  PUT  , OUTPUT  ,TApF3  ) 

CO‘^HriJ/ITl/IT,  TOPIC,,  ncLj  ,Tnu-%X'(  ‘^152) 

OTHCrSTCM  ISUH(<iC)  , J^U'' ( 6?  ) , XT”  ( ?1 .3 ) 

EQUlVfLEMOE  (XTT(1),TT)  , , ( HCLT  , 0"LTI‘M 

PEfiO^  ,CLKPT  ,SaMPPC  jCARP^n,  MSA^P,  U,  I^SK  ,NSFO  , TORIG 

MTTL=  ( -1 

nELTIH=l/SAMFRO 

nEL?CFT=l/CLKRT 

PI2  = r . 28312  530  8 

TIKF=C  .0 

SHIF->-  = C.  0 

no  1:0  I = 1,HPA'''° 

Tr(TIMF.GF.RHirT)GO  TO  101 
GO  TO  102 

101  CALL  P^(^^‘!AS<,MG^p,^'P:T,NCNT) 
ISUM(Nr!,|T)=IFU'^(|.'CNT)  +1 
GHIF‘"=r.HIFT  + DELEHPT 

1 C2  VT  ( I)  = ( (Nr:NT-(N/?)  )*  (COS  (PT2^CARFR')»TI'^F)  ) ) / (*|/2) 

100  TI'IFrTTME  + OfrLTIH 
L'N=rjr  APP  + 1 
WRITP(3)  XTT 

IF  ( I'".  GT.  3 12)  WRITE  (3)  ( XT  ( J ) , J=  "1 3,  T T ) 

WRITE  1,IT 
VIRITP  2 

WRITP  3, (I, ISUM(I)  ,I  = 1,N) 

WPITP  u 
STOP 

1 FOPwrTf'^  DATA  STORFO  O'J  OATA  SET  3 tJPT?  = ‘,IR) 

2 FOP‘’(T(’  STATISTICS  OP  0ATf*,./,'‘  WCNT^) 

FOPMfT  (0v(lX,I?,3X,TP,/)  ) 
cQOMrTC*  STATISTICS  OF  PON*,/,«‘  NCNT») 

5 FOR^^f  T (GO  (IH  , 12 , 3 X , I 8 , / ) ) 

FNO 


Figure  38.  SPRNGEN  Main  Program 
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lOFNT  3VJ 


jODDfii  — 


1R(Mt 


SUr^^^OlJTT'JF  TO  GrNEf?tTr  a os?:Lno  NOIG-  (P'i) 
CALL  PMcrulASKjMS-Ot'niT,  UCfJT) 


M 

LENGTH  OF 

TME  SE0'Jc-;jcE  IN  BITS. 

MASK 

'^-FORA  CK 

K'S'^n 

SEOUFT'CE,  ‘^'JST  RE  INITIALI7E0  N 

NOTT 

RI‘-HT>-'05T 

PIT  0-  N'REP. 

NC‘'T 

NJHRER 

ONES  in  NSEO. 

ENTRY 

RN 

RA'^F 

0 

USE 

CDO'" 

A1 

CONTAINS 

AROPESS  OF  APGU’^ENT  LIST. 

'^SS 

1 

SA? 

A 14-1 

r,PT  MA$'<  IN  X2. 

SA? 

X2 

SA7 

A14-2 

GET  NSEO  IN  X3. 

SA3 

X3 

BX2 

X2'‘X3 

logical  Atn. 

CX2 

X2 

CO'JNT  ONES  IN  XR  . 

mxF 

-1 

GENERATE  A ONf  jn  XG. 

RX*" 

-XE 

RX? 

X?‘f-XS 

EV-N  OR  nno  NUMBER  0'=’ 

St^h 

XI 

GET  N IN  XL. 

SRi- 

VJ*. 

LX2 

R'^ 

SHIFT  LE*^!  N PITS. 

BX"’ 

X3+X? 

INCLUSIVE  OP. 

AX7 

1 

shift  right  OMr  on. 

SA7 

A3 

R'^THRN  NEW  NSEO. 

PX«-- 

XS*^^ 

LOGICAL  ANO, 

SA? 

A 1 + 3 

RETURN  NEW  NRIT. 

SAF 

X2 

CX7 

X7 

NUNPPR  OF  ON>^W  IN  NFEQ 

SA? 

A 1 + •* 

RETURN  NEW  NCNT, 

SAR 

X2 

JP 

ENP 

•RN 

RFTflRN. 

Figure  39-  Assembly  Language  Subroutine 
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